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What are the fundamental laws for the adsorption of charged polymers onto oppositely charged 
surfaces, for convex, planar, and concave geometries? This question is at the heart of surface coating 
applications, various complex formation phenomena, as well as in the context of cellular and viral 
biophysics. It has been a long-standing challenge in theoretical polymer physics; for realistic systems 
the quantitative understanding is however often achievable only by computer simulations. In this 
study, we present the findings of such extensive Monte-Carlo in silico experiments for polymer- 
surface adsorption in confined domains. We study the inverted critical adsorption of finite-length 
polyelectrolytes in three fundamental geometries: planar slit, cylindrical pore, and spherical cavity. 
The scaling relations extracted from simulations for the critical surface charge density a c —defining 
the adsorption-desorption transition—are in excellent agreement with our analytical calculations 
based on the ground-state analysis of the Edwards equation. In particular, we confirm the magnitude 
and scaling of a c for the concave interfaces versus the Debye screening length 1/k and the extent 
of confinement a for these three interfaces for small na values. For large na the critical adsorption 
condition approaches the planar limit. The transition between the two regimes takes place when the 
radius of surface curvature or half of the slit thickness a is of the order of 1/k. We also rationalize 
how <t c (k;) gets modified for semi-flexible versus flexible chains under external confinement. We 
examine the implications of the chain length onto critical adsorption—the effect often hard to tackle 
theoretically—putting an emphasis on polymers inside attractive spherical cavities. The applications 
of our findings to some biological systems are discussed, for instance the adsorption of nucleic acids 
onto the inner surfaces of cylindrical and spherical viral capsids. 


I. INTRODUCTION 


The adsorption of charged polymers orpolyelectrolytes 
(PEs) onto oppositely charged surfaces has a num¬ 
ber of technological and biophysical applications includ¬ 
ing paper production & H] , interface coating |[10|, layer- 
by-layer formation [3 7], water desalination [ill Il2| , and 
stabilization of colloidal suspensions 0 , Il4j. One dis¬ 
tinguishes weak and strong PE-surface adsorption 
weak PE-surface adsorption is governed by an interplay 
of energetic often electrostatic (ES) attraction of polyions 
onto an interface versus an entropic penalty accompany¬ 
ing the confinement 0,113 or the compression of the 
polymer chains near the surface 0 - 0 ]. Weak adsorp¬ 
tion takes place for weakly charged PEs (partially neut¬ 
ralized by condensed counterions [2tl IzI ) and for mod¬ 
erately charged interfaces: such PE-surface binding is 
rather reversible. This contrasts an irreversible adsorp¬ 
tion in the limit of strong PE-surface association [g, |23| . 


Weak PE adsorption onto surfaces of different geomet 
ries at varying conditions has been investigated in a num¬ 
ber of recent theoretical 15 


and computer simulation 


2J-3 


52T6S 


j, experimental j40H51j . 
studies. PE adsorption 


onto patterned and corrugated surfaces [30J, [3l|, |70H73( | as 
well as critical PE adsorption onto charged Janus net- 
neutral particles (52 ] was also examined. The proper¬ 
ties of polymer adsorption inside cylindrical nanopores 
in porous glass were studied experimentally and by com¬ 
puter simulations in Ref. 0] . Possible effects of surface 



Figure 1: Schematics of the inverted PE adsorption in confine¬ 
ment: planar slit, cylindrical pore, and spherical cavity. Video 
files illustrating the conformational changes of the polymer 
chain for conditions below and above the critical adsorption 
transition are presented in the Supplementary Material. 


dielectric discontinuities on PE-surface adsorption were 
rationalized too 0-0 . 

The critical adsorption describes the threshold condi¬ 
tions at which the ES-driven adsorption of PE chains 
first takes place in the system. This phase transition in¬ 
terrelates the condition for the interface surface charge 
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density er, the line charge density of the polymer p, the 
reciprocal Debye screening length in the solution k, the 
ambient temperature T, and the polymer’s Kuhn length 
b. The critical adsorption condition defines the relation 
between these important model parameters at the co¬ 
existence boundary of adsorbed versus desorbed chain 
conformations. Typically, a universal critical adsorption 
parameter can be constructed, 

_ 247ra 3 |pcr e | 

c ek B Tb ’ 1 ’ 

and its dependence on k governs the scaling of cr c —the 
critical surface charge density required for PE adsorption 
to take place—at varying salt conditions. Here e is the 
dielectric constant of the medium, a is the curvature ra¬ 
dius of the adsorbing surface, and ks is the Boltzmann 
constant. For a long flexible nearly Gaussian charged 
polymers in front of a uniformly oppositely charged plane 
the well-known result is E2 

#(k)~^(k)~« 3 - (2) 

The standard rationale for increasing a c with the salt 
concentration no is the requirement to compensate a 
stronger screening of ES attraction of the PE chain to 
the oppositely charged surface. For a symmetric 1:1 elec¬ 
trolyte we have n 2 = 8nlsno : where Ib = ^/(eksT) is 
the Bjerrum length. The peculiar cubical scaling of the 
critical charge density with n in Eq. m stems from the 
properties of the eigenfunctions of the corresponding Ed¬ 
wards equation for the conformations of a long polymer 
chain in the attractive Debye-Hiickel potential of the in¬ 
terface E3- In addition, some ES chain stiffening at 
low-salt conditions takes place impeding the PE-surface 
adsorption (see also Ref. [H] for non-ES effects in PE- 
surface adsorption). 

For the convex cylindrical geometry (see Fig. [TJ a 
quadratic scaling is instead predicted by the Wentzel- 
Kramers-Brillouin (WKB) theory at low salt [3(|, namely 

Sf(na) ~ M 2 , (3) 

while at high salinities and large rod radii m > 1 the 
planar limit © is recovered. For PE adsorption on the 
outside of oppositely charged spherical particles yields 
the linear dependence of cr c («;a) in the limit 1 [36}, 


micelles of spherical and cylindrical geometry [43| . These 
experimental findings indicate a weaker dependence of 
a c on na for more ’’convex” surfaces, as the adsorbing 
interfaces transfer from the planar to the cylindrical and 
finally to the spherical shape. 

The adsorption transition of weak PEs under con¬ 
finement [29[—we call below the inverted critical 
adsorption— has a number of biologically relevant applic¬ 
ations. Here, the term critical for confined PEs has the 
same meaning as for adsorption of charged polymers onto 
the planar and convex interfaces For instance, the 
self-assembly of cylindrical and spherical single-stranded 
RNA viruses involves the adsorption of nucleic acids onto 
the inner virus capsid surface composed of protein build¬ 
ing blocks [23, [33^ LZ3 t [Z^|. The capsid proteins are abund¬ 
ant in liighly-basic flexible poly-peptide tails which trig¬ 
ger the adsorption of negatively-charged nucleic acids 
thus steering the self-assembly [73, [8(| • The known ex¬ 
amples include nucleic acid encapsulation inside the cyl¬ 
indrical tobacco mosaic virus TMV [8l| an d the icosa- 
hedral cowpea chlorotic mottle virus CCMV [27l l82[ i83}. 

For very long chains, the scaling relations for critical 
PE adsorption in inverted geometries were recently de¬ 
rived theoretically from the ground-state analysis of the 
Edwards equation for the Green function [37j . The main 
subject of the current paper is the inverted weak adsorp¬ 
tion of finite-length PE chains of varying stiffness in all 
three basic geometries, see Fig. [Q For long flexible poly¬ 
mers the critical adsorption conditions were obtained us¬ 
ing the WKB method in Ref. j36|. Namely, Eqs. (5), 
(6), and (7) of Ref. [36] provide the dependence of a c in 
the entire range of na. At low salt or strong confinement 
when na < 1 the critical adsorption parameter 8 l ™ for 
a planar slit was predicted to scale as [37] 

^pl,inv ~ 3C '2( Ka )l. ( 5 ) 


Here and below the constant C is of order unity, C ss 
0.973. For long flexible PEs inside the oppositely charged 
cylinder this parameter reveals a plateau with a slowly 
varying logarithmic correction [ 37 J 


S cyl 


,mv 


3 c 2 

0.116 — log(«;a) 


( 6 ) 


Finally, for a PE inside spherical cavities in the low salt 
limit the value of 5 C tends to saturate to a plateau [37} 


5 s c p (na ) ~ (no) 1 . 


( 4 ) 


^sp,inv ^ 0^2 


( 7 ) 


Here a is the radius of the cylinder or sphere. For more 
details on these scalings the reader is referred to Refs. 


The systematic change in the <r c (rea)-scaling behavior 
from the planar interface via a cylinder to a sphere is in 
agreement with a number of experimental evidences from 
the Dubin’s lab, see e.g. Ref. [43| (and also the analysis 
in Ref. HI)- The experimental observations of critical 
PE adsorption are based on the complex formation of 
various polymers with oppositely charged particles and 


These functional dependencies on na are in stark con¬ 
trast to the fast and monotonically increasing S c for the 
adsorption of PEs on the outside of cylindrical and spher¬ 
ical interfaces, Eqs. © and ©. In the opposite limit of 
loose confinement or high salt when na 1 the theory 
predicts 

§ inv = 3 C 2 (ko ) 3 /2 (8) 

for all three inverted geometries [37]. This latter limit¬ 
ing behavior was derived from the WKB approach in the 
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limit of zero surface curvature H3- Note that for a 
finite-length polymer all the above mentioned standard 
ground-state based predictions for a c need to be mod¬ 
ified; the regular procedure is however not easy. This 
makes the findings of our computer simulations—the 
main focus of the current investigation—even more valu¬ 
able for experimentally relevant situations. 

In the current paper we study by extensive Monte- 
Carlo simulations the properties of critical PE adsorption 
in three basic inverted geometries, see Fig. [T] We study 
the effects of the chain length, the polymer persistence, 
and systematically of the confinement size and solution 
salinity onto the critical surface charge density a c . In 
Sec. m we present the details of the simulation model 
and the data analysis algorithms. The main results on 
PE adsorption profiles and critical adsorption character¬ 
istics are described in Sec. uni We discuss the physical 
rationales behind the observed dependencies and the ap¬ 
plications of our results in Sec. ED 


II. MODEL AND APPROXIMATIONS 


We here implement the same Metropolis Monte-Carlo 
simulation algorithm which has been successfully ap¬ 
plied and tested by us recently for PE adsorption onto 
spherical (H3| , cylindrical [HU, and spherical Janus [H2] 
particles as well as for surface adsorption of pH-sensitive 
PEs IfiH. We refer the reader to Refs. |52l Idfl . [57} 
for more details on the simulation procedure. In brief, 
the polymer chain is modeled within the spring-bead 
model, with each monomer being a rigid sphere of radius 
R m = 2A carrying a point elementary charge Zm — 
at its center. Neighboring beads are connected by the 
harmonic potential 

U str (r ) = K r (r - r 0 ) 2 /2, 

with the elastic constant for bond stretching K r = 
l.ON/rn 2 and the inter-monomer equilibrium distance 
ro = 7.0A (as for single-stranded DNA 0). The chain 
stiffness is given by the elastic potential 

U ei {9) = K e (6 - 9 0 ) 2 /2, 


where the force constant Kg assumes the values such that 
the non-ES persistence length l p $ of the polymer ranges 
from about 8 to 50 A (a typical range for many real 
PEs). Here 6 denotes the angle between the two success¬ 
ive bonds and 9q = n. The mechanical persistence length 
for an uncharged chain l p ^ was obtained in simulations 
via the relation l Pt o = (i? 2 ) 1 / 2 /(l + (cos0)), where (R 2 ) is 
the root-mean-squared monomer-monomer distance [85|. 
The inter-chain excluded volume is accounted for by the 
standard hard-core repulsive potentials in simulations, as 
compared to the theoretical model [13]. 

The repulsion of monomers at distance r is given by 
the screened Coulomb potential 


Ues{t) 


er 


(9) 




Figure 2: (A) Variation of the dimensionless ES potential 
'l'(a) = eo ifi(a)/(ksT) at the boundary versus the radius 
of the spherical cavity plotted for: a = — 0.1/(4-7r) C/m 2 
and k=1/(30A) (filled circles), a = —0 .1/(87 t) C/m 2 and 
k=1/(30A) (filled squares), a = —0 .1/(167r) C/m 2 and 
k=1/(30A) (filled triangles), a = —0 .1/(167r) C/m 2 and 
k= 1/(100A) (empty triangles). Note that C/m 2 « eo/(16A 2 ) 
that is ~ 6 <tb-dna, where ctb-dnaIs the bare charge density of 
the B-DNA [84;]. For large a values the potential approaches 
the known planar result, ^p(a) = 47 r<r/(eK). (B) The ES poten¬ 
tial on the inner surface of spherical cavities scales as 4/ (a) ~ 
na for na 1, whereas ^(a) —> eo<fi{a)/(I cbT) for Ka 1. 
Parameters: a = — 0.1/(4-7r) C/m 2 , k = 1/(30(X)A) (blue sym¬ 
bols), k = 1/(1000A) (red symbols), k = 1/(300A) (green 
symbols). 


The ES potential emerging in a slit with inter-plane 
distance 2a, inside a cylinder or a sphere of radius a 
were computed as the solutions of the linear Poisson- 
Boltzmann equation [37}. Below, we use the potentials 
denoted as Tm,oiit( r ) i n [37} to parametrize the 

strength of ES PE-surface attraction. For brevity, we 
do not provide the explicit analytical expressions here, 
instead showing the potential distributions in Fig. [2] 
The critical surface charge density a c is defined in our 
simulations as the condition at which the PE binding 
energy to the interface exceeds the thermal energy, 

\E b \ > k B T. (10) 

Thus, even for the conditions when the polymer is not in 
direct contact with the surface but its total binding en¬ 
ergy is lower than E b = —k B T , we consider the chain to 
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be in the adsorbed state. To compute the value of a c for 
given values of the model parameters k, a , and N, we per¬ 
form the simulations for a set of surface charge densities 
and then determine the one for which the adsorption- 
desorption criterion m is satisfied. Here one can an¬ 
ticipate already that longer chains will require smaller 
surface charge densities cr c to be classified as adsorbed, 
as we indeed obtain from simulations, see below and also 
Ref. [86|. In the limit a —> oo the potentials ’P iraiOTit (r) 
of Ref. [33 turn into (r) for the corresponding isolated 
surfaces (see Fig. ED; the same holds for the properties 
of the adsorption-desorption transition, see below. 

There exists a number of differences between the in¬ 
verted PE adsorption and the polymer-surface adsorp¬ 
tion from a dilute, free-space solution. One feature is 
the presence of confining interfaces. They have differ¬ 
ent implications onto the polymer: for a planar slit the 
polymer is mobile in two dimensions, for a cylindrical 
tube the chain is free to move in one direction, and for 
a spherical cavity the polymer has no translational free¬ 
dom at all. This progressively increasing confinement 
reduces the polymer conformational entropy [871 ] . par¬ 
ticularly upon adsorption on the interior of oppositely 
charged cylinders and spheres, see also Ref. |8a |. 

We also note that in the low-salt limit the total PE 
persistence length, 

l P (K) = l P ,o + l* S (K) = b/ 2, (11) 

acquires an ES component which is decreasing with the 
solution salinity. For flexible chains it obeys the scal¬ 
ing 1% s (k) ~ K - 1 i while for semi-flexible polymers 
~ K- 2 Q. This fact is not accounted for in the 
theories of PE-surface adsorption [H,[37j yielding for 5 l J lv 
the scaling relations @, @, 0. This ES contribution 
lp S should renormalize the scaling of cr c with k obtained 
from computer simulations at low salt, in accord with 
Eq. 0, see Ref. [HI the for planar and convex surfaces. 
For concave adsorbing interfaces, such as the sphere’s in¬ 
ner surface, due to this ES polymer stiffening the chains 
will tend to occupy regions of smaller curvature [8?| . as 
it is indeed observed upon ’’spooling” of double-stranded 
DNA inside bacteriophages [89(. The chains approach the 
interface because of bending energy minimization thus 
facilitating the ES-driven PE-surface binding, see below. 

An additional important parameter for confined PE- 
surface adsorption is the polymer’s volume density. In 
the theory y|, |3f| the PE adsorption typically takes place 
from a very dilute polymer solution, which is not the case 
for confined inverted-adsorption situations, where the net 
polymer density is finite, see below. 

Let us now briefly discuss some approximations in¬ 
volved in the current study. 

a) We use the Debye-Hiickel theory to compute the ES 
potentials near the interfaces and between the polymer 
monomers. This approach is valid for weakly charged 
systems and for an appreciable amount of salt in the 
solution no, when the ES potentials |'F| <25 mV (com¬ 
pare the panels of Fig. 0. The solution of the nonlinear 


Poisson-Boltzmann equation in curved geometries in the 
presence of salt is a formidable theoretical problem per 
se, often only solvable in some idealized limits (see Ref. 
[901 ] for a charged rod at n 0 = 0). Note that the linear 
ES theory often overestimates the magnitude of the po¬ 
tential emerging near highly charged interfaces (see Fig. 
2 of Ref. (91| and Ref. [S3). Also note that partic¬ 
ular in low-salt solutions, the effects of counterion re¬ 
lease from the surface—on the level beyond the standard 
Poisson-Boltzmann approach with the cation concentra¬ 
tion obeying n(r ) = n^e —might become relevant 
for PE-surface adsorption. 

b) The WKB scaling relations presented in the Intro¬ 
duction stem from the ground-state analysis of infinitely 
long flexible Gaussian chains in front of surfaces with 
the Debye-Hiickel ES potential. Both these idealizations 
will not hold upon variation of no i n a broad range, as 
we study below. The investigation of implications of the 
non-linear nature of the ES potential near highly charged 
surfaces is the subject of a separate investigation [931 ]. 
Also, the impact of the mutual influence of adsorbing 
PE chains onto the ES potential of the interface (charge 
regulation) can non-trivially impact the critical adsorp¬ 
tion conditions in terms of <j c (k) scaling. Moreover, we 
consider below the adsorption of a single PE chain from 
a dilute solution; at realistic conditions however several 
chains might adsorb simultaneously. Their mutual salt- 
dependent ES repulsion along the surface will have an 
effect, for instance, on the overall surface coverage by 
PEs. The latter is often measured experimentally for 
PE-surface adsorption from bulk solutions with a finite 
polymer concentration. All these effects are experiment¬ 
ally relevant and will be considered elsewhere [93| . 

c) We implement the /csT-based adsorption criterion 
(fTUl) to identify the PE adsorption-desorption threshold. 
This criterion is however somewhat arbitrary: for in¬ 
stance, one can classify the adsorption threshold us¬ 
ing polymer distributions with a single peak versus the 
double-peaked profiles emerging between the two confin¬ 
ing interfaces, see e.g. Fig. [3]below. We assume reversib¬ 
ility and ergodicity [94| for the process of PE adsorption 
at all conditions (no irreversible binding). This assump¬ 
tion might not be valid, particularly at low salt when 
the binding of even several PE monomers to an interface 
with large ES potential might overcome the thermal en¬ 
ergy. Also, in this limit the relative accuracy in defining 
the adsorption-desorption boundary becomes important 
and even small fluctuations |<5ct c |/|<t c | that can realize 
in experiments might cause sizable effects. Finally, the 
adsorption of one fragment of the chain is assumed not 
to affect the ES potential acting to attract other parts 
of the polymer. This might be important for pH- and 
potential-responsive surfaces (not a part of this study 
[93|). We work in the single-chain limit and thus do not 
study the PE adsorption isotherm—the amount of poly¬ 
mer adsorbed for a varying bulk polymer concentration. 

Despite these simplifications and assumptions, our 
computational results reveal excellent agreement with 
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Figure 3: Distribution of polymer monomers for inverted ad¬ 
sorption of flexible PEs (l p ,o = 8 A.) in a spherical cavity (top 
panels), cylindrical pore (middle panels) and planar slit (bot¬ 
tom panels) with the surface charge density of a = — 0.1/(47t) 
C/m 2 . The radius of the sphere and cylinder is a = 50 A, and 
the slit thickness is 2 a = 100 A. The distance r denotes the 
separation from the centre of the confining space. The degree 
of chain polymerization is IV = 20 (left panels) and N = 100 
(right panels). The salt concentration no is varied: kb = 0.1 
(black), na = 0.5 (red symbols), na = 1 (blue symbols) and 
K,a = 5 (green symbols). 


the theoretical predictions in a wide range of model para¬ 
meters, see below. 


III. RESULTS 

A. Polymer density distribution 

First, we examine the distribution of polymer 
monomers, p(r), in the three basic inverted geometries. 
For a fixed degree of the external confinement, the evol¬ 
ution of the p(r) profiles with varying no reflects the 


Figure 4: Distribution of semi-flexible polymers inside a 
spherical cavity (top), cylindrical pore (middle) and planar 
slit (bottom panels) for the parameters of Fig. [3] and for 
N = 100, na = 0.5 (left panels, reddish colors) and na = 5.0 
(right panels, greenish colors). The color scheme matches that 
of Fig. [I!] The non-ES persistent length is l Pl o = 8 (circles), 
23 (squares), 33 (diamonds) and 46A (triangles). 


positioning of the adsorption-desorption boundary. We 
start with relatively flexible chains confined into a spher¬ 
ical cavity and a cylindrical pore; we use l Pi o = 8A for 
flexible chain results below. If the unperturbed radius 
of gyration of the chain exceeds the cavity dimensions, 

(Rg) > a i volume exclusion creates a force pushing 
the polymer towards the surface, see Fig. [3] This is par¬ 
ticularly pronounced for spherical cavities, with strongly 
restricted chains. With increasing ionic strength up to 
na — 1 the monomer accumulation near the surface gets 
facilitated. From na = 1 to na = 5 this behavior gets 
inverted, as at large na the conditions are close to or 
above the threshold of polymer desorption, and in ad¬ 
dition the ES term in the chain persistence length gets 
smaller. This makes the polymer chains effectively more 
flexible, they are attracted weaker to the interface, and 
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the polymers occupy the bulk of the cavity more readily. 
In this situation we thus find a single peak of the polymer 
distribution in the centre of the confining space, whereas 
for small na the PE peak emerges near the oppositely 
charged surface in confined geometries. 

The effect of the polymer’s mechanical persistence onto 
PE distribution in confined spaces is illustrated in Fig. [4] 
Progressively stiffer PE chains prefer to occupy the peri¬ 
pheral regions of the cavity due to a lower bending energy 
penalty for the chain arrangements with larger radii of 
curvature - the effect particularly pronounced for spher¬ 
ical cavities, see the top row in Fig. [I] These latter trends 
are in line with the results of our recent simulations for 
polymer chains inside inert or non-attractive spherical 
cavities in the presence of macromolecular crowding [ 83 ] , 
see also Refs. [H, [9(|. Moreover, the loss of configur¬ 
ational entropy for the arrangement of more persistent 
chains near the cavity surface is smaller as compared to 
the flexible ones. Fig. [4] shows that the deviations for 
semi-flexible versus flexible chains become smaller as we 
go from PE adsorption inside a spherical cavity to PE 
adsorption inside a cylindrical pore and finally to PE ad¬ 
sorption inside a planar slit (respectively, the top, middle 
and bottom panels of Fig. 11. The physical reason is 
again the the number of the polymer’s degrees of freedom 
available in the corresponding geometries. For instance, 
for a planar slit the bending energy of semi-flexible chains 
has nearly no implications on the amount of the poly¬ 
mer near the adsorbing interface, in stark contrast to 
the spherical cavity for which a severe chain bending 
is unavoidable [ 13 ], compare the panels in Fig. d As 
we demonstrate below, this polymer bending energy in 
spherical and cylindrical confinement has non-trivial ef¬ 
fects onto the critical surface charge density a c (Ka) of the 
adsorption-desorption transition. 

A relevant experimental question for super-critical PE 
adsorption is the amount of the polymer adsorbed on the 
surface. In Fig. [5] we quantify how this amount changes 
with Ka for the three adsorption geometries for the single¬ 
chain adsorption simulated. Namely, for a > a c we ana¬ 
lyzed the PE profiles formed near the interfaces for the 
inverted polymer adsorption. We evaluate the fraction 
of the polymer chain in the region close to the adsorbing 
interface, N a d s /N. This fraction is a stationary quant¬ 
ity; we do not consider here the kinetics of PE adsorp¬ 
tion for a > a c H3, HI] ■ We illustrate the behavior 
of this fraction in Fig. [5] versus the reciprocal Debye 
screening length, n. We find that the amount of PE 
adsorbed within this first layer near the interface is of¬ 
ten a non-monotonic function of the salt concentration 
n o- One physical reason for this is a shorter ES persist¬ 
ence length of PEs and weaker polymer-surface charge- 
mediated binding as no increases. This non-monotonicity 
can be anticipated already from the evolution of PE pro¬ 
files in Fig. [3] in the proximity of the adsorbing sur¬ 
face as k increases. Here, we refer the reader to the 
studies in Refs. [13, E| for experimental evidence and 
theoretical predictions of non-monotonic effects of ad- 



Figure 5: Fraction of polymer monomers within 4A from the 
adsorbing surface for inverted PE adsorption into a spher¬ 
ical cavity (black), cylindrical pore (red) and planar slit (blue 
symbols) with a = — 0.1/(47r)C/m 2 (above the critical ad¬ 
sorption transition), plotted for varying solution salinity. The 
radius of spherical cavity and the cylindrical tube is a = 50 
A and the inter-plane distance for the slit is 2 a = 100 A; 
N = 20 (left panel) and N = 100 (right panel); l P: o = 8A. 


ded salt on the amount of adsorbed PE chains. In real¬ 
istic multi-chain systems, higher solution salinities effect 
softer PE chains, weaker PE-surface ES attraction, but 
also a weaker ES repulsion between the already adsorbed 
polymer coils. The interplay of these effects might yield a 
non-monotonic behavior of the mass-per-area of adsorbed 
PEs with varying k. 

The width of the PE profile w together with the mass 
of adsorbed PE per area are the experimentally relev¬ 
ant quantities. As we work in the single chain limit the 
latter will not be considered. For inverted PE adsorp¬ 
tion, the width of the adsorbed layer is expected to be 
a non-monotonic function of k too. Note however that 
the standard definition of w implemented for instance for 
the adsorption of an isolated PE chain onto an attractive 
surface [5[—as the width of the polymer probability dis¬ 
tribution function at its half-height—cannot be directly 
used for the current problem of inverted adsorption. 


B. Critical adsorption conditions 

Now we turn to the main objective of the current study, 
the scaling of the critical surface charge density at the 
adsorption-desorption transition. For the inverted PE 
adsorption in confined geometries, we find that a c var¬ 
ies with Ka as shown in Fig. [6l revealing an excellent 
agreement with the theory developed in Ref. j]37] . As ex¬ 
pected, the scaling of the critical surface charge density 
is very different in the limit of small and large Ka values, 
as prescribed by Eqs. ©, ©, © and Eq. ©, respect¬ 
ively. In the low-salt limit spherical cavities necessitate 
larger surface charge densities than the cylindrical tubes; 
the latter in turn need larger surface charges than the 
planar slits in order to reach the same degree of polymer 
binding (fTOl) . compare the curves in Fig. [Gj We attribute 
this reduction of a to a progressively smaller penalty 
of entropic confinement of flexible PE chains. Remark- 
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Figure 6: Critical surface charge density a c for inverted ad¬ 
sorption of flexible PEs (l Py o = 8A) into a planar slit (blue 
symbols), cylindrical tube (red symbols), and spherical cav¬ 
ity (black symbols), plotted for varying solution salinity. The 
dashed lines are the full theoretical asymptotes for 5 C [37i |. 
plotted for a c using Eq. m- The corresponding low-salt 
limits are given by Eqs. ©, ©, 0 ; the high-salt or the 
planar limit is Eq. 0. Note the inverse positioning of the 
curves in the region of /to S> 1. Parameters: a = 50 A (1/2 of 
the slit thickness, the cylinder and sphere radii), the polymer¬ 
ization degree is N = 20 (panel A) and N = 100 (panel B). 
On a standard 3-3.5 GHz workstation every curve on these 
graphs requires some 180 h and 900 h of computational time 
for chains of N = 20 and N = 100 monomers, respectively. 


ably, even the reversed order of the critical adsorption 
curves in the limit of large /ta, as compared to the low- 
salt limit, is precisely reproduced in our simulations, in 
accord with the theory [371 ]. This reversed order at high 
salt is particularly well pronounced for longer chains, Fig. 
[6)3. The change in the cr c (na) scaling behavior, from the 
low-salt prediction to the high-salt asymptote (0, oc¬ 
curs at na ~ 1, for any chain length, compare the panels 
in Fig. [6] The effects of the surface curvature on the 
adsorption-desorption transition point is thus universal 
and start to be important for the conditions of low-salt 
and large surface curvature when n < 1/a. Fig. [6] is the 
central result of the current study. 

For progressively longer PE chains the values of a c de¬ 
crease, preserving however the overall scaling relations 



Figure 7: The same as in Fig. [0 for inverted adsorption inside 
a spherical cavity (black), cylindrical tube (red), and planar 
slit (blue symbols) for varying re. The results are plotted in 
the linear scale. The polymer stiffness is l Pt o = 8 A (open 
symbols) and l Pt o = 50A (full symbols). Parameters: a = 54 
A and N = 20. 


in the limit of low salt and the approach to the univer¬ 
sal asymptote at high salt, compare Figs. 0A and 03. 
According to our adsorption criterion, Eq. m, longer 
chains will clearly accumulate the same binding energy 
for smaller |cr|. A similar behavior was observed for PE 
adsorption onto Janus particles [52|: the entire a c (na) 
dependence shifts down for longer chains but preserves 
the scaling properties for small and large rea values. 

In Fig. 0 flexible chains are considered, with the non- 
ES persistence length of l Pi o < 10 A. For a spherical 
cavity, the theoretical results of Ref. !37j for a c give an 
excellent agreement with the results of our simulations 
for N = 100 chains if we set b = 2Z p .o «4 A, a realistic 
value for flexible chains simulated. This follows from Eq. 
0 for PEs inside spherical cavities for na <C 1 from 
which the critical charge density can be recalculated in 
units C/m 2 as 


sp,inv 

° c _ 

C/m 2 


2 16fer 0 
247ra 3 Z B 


1.9 x 1CT 5 . 


( 12 ) 


The same Kuhn length b is then used to compute the full 
a z c nv (n) asymptotes from the theoretically predicted [37} 
relations for 8^ yl ’ mv (n) and 8 yl,mv (n), according to Eq. 

0 - 

Due to the adsorption criterion implemented, for N = 
20 chains a c is nearly 100/20=5 times larger than for 
N = 100 polymers, compare the panels in Fig. [6] and 
see also the universal curves in Fig. [lU] Also note that 
for na 1 our simulations of the PE adsorption under 
confinement yield the a c (n) ~ re 3 scaling behavior, as an¬ 
ticipated for polymers with a salinity-independent Kuhn 
length b E3- This is in contrast for instance to PE ad¬ 
sorption on the outside of spherical particles, where the 
effects of ES persistence are important and our simu¬ 
lations in the high-salt limi t gi ve cr c (n) ~ re 1 ’ 9 scaling 
instead, see Fig. 8A in Ref. [52 ]. 

We also examined the dependence of the critical ad¬ 
sorption conditions for more persistent chains in all three 























Figure 8: Critical adsorption charge density for persistent 
versus flexible polymers inside a spherical cavity, cr c /cr[?, com¬ 
puted for varying non-ES polymer persistence length. Para¬ 
meters: a = 50 A, N = 100, K,a = 1 (circles), 5 (squares), 8 
(diamonds) and 10 (triangles). 



Figure 9: The same as in Fig. [6] but for varying radius a of 
the spherical cavities. The high-salt asymptote © for the 
charge density (in units of C/m 2 ) is obtained from Eq. © 
as cr c (K)/[C/m 2 ] = C 2 K 3 bro/(nlB) ~ k 3 , the dashed line. 
Parameters: N = 100, a = 50A, l p , o = 8A. 


adsorption geometries, see Figs. [7] and [8] Fig. [7] shows 
that for more persistent chains the magnitude of a c de¬ 
creases for the adsorption onto the planar slit, inside cyl¬ 
indrical pore, and spherical cavities. We find that par¬ 
ticularly for PE adsorption inside spherical cavities the 
magnitude |cr c | decreases due to a bending-energy driven 
localization of polymers near the cavity surface, see Fig. 
[5] For PE-sphere inverted adsorption, stiffer PEs prefer 
to stay closer to the adsorbing interface thus reducing 
the value of a c . The precise behavior of u c as a function 
of k a shows that the deviations from the flexible chain 
results become progressively larger for more persistent 
chains and at larger na values, see Fig. © The latter 
is not surprising because at high-salt conditions the ES 
contribution to the polymer persistence gets reduced and 
the PE stiffness is dominated by its mechanical part l p ^, 
see Eq. (fill) . 



Figure 10: Universal rescaled surface charge density for crit¬ 
ical PE adsorption inside spherical cavities of varying radii 
and for polymers of different lengths, plotted for varying k; 
l P , o = 8A. 


The question arises whether for inverted PE critical 
adsorption the variation of the confinement degree or 
salt concentration gives rise to different a c behaviors? 
We showed that for adsorption of finite-length PEs onto 
spherical Janus particles [52[ there was no universal para¬ 
meter na that would combine the curvature and salinity 
effects on the critical adsorption properties cr 0 (/ta). The 
inverted critical PE adsorption is also quite different if 
one varies the size of the confined cavities, tubes, and 
slits or the solution salinity, compare the curves in Fig. [9] 
for spherical confinement. Longer chains require smaller 
surface charge densities to get adsorbed and for smaller 
sizes of spherical cavities the value of a c increases. This 
can be understood from the variation of the ES surface 
potential in spherical cavities presented in Fig. ©a show¬ 
ing that iff (a) « 2(Ka)eoip(a)/(ksT) ~ na for na <C 1 
|37| . For inverted critical PE adsorption we obtain that 
indeed there exists no universal parameter na. This at 
first sight disagrees with the theoretical results of Ref. 
Ifl'fl . For the finite-length PEs with varying ES persist¬ 
ence studied in our simulations this disagreement is how¬ 
ever not surprising, as compared to infinitely long flexible 
salt-insensitive polymers studied in the theoretical ideal¬ 
ization [ 1133 - 

Fig. [Tell illustrates the behavior of the rescaled critical 
surface charge density for inverted PE-sphere adsorption, 
a c aN. This combination accounts for the peculiar fea¬ 
tures of the variation of the ES potential with the cavity 
radius (Fig. ©r) and the adsorption condition used in 
the simulations, Eq. m■ We find a universal collapse 
of this renormalized parameter for the critical adsorption 
curves for different chain lengths and cavity sizes, for a 
wide range of variation of the solution salinity. 


IV. DISCUSSION 

In this study, we have employed extensive Monte-Carlo 
computer simulations to unveil the physical properties 
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of PE adsorption in confined spaces, considering poly¬ 
mer chains inside a planar slit, a cylindrical pore, and 
a spherical cavity. We rationalized the position of the 
adsorption-desorption transition upon variation of vari¬ 
ous physical parameters such as the extent of the ex¬ 
ternal confinement, the salinity of the solution, the chain 
length, and the bare persistence length of the polymer. 
We have demonstrated how the well-known cubic scal¬ 
ing of the critical surface charge density with the recip¬ 
rocal Debye screening length n gets non-trivially mod¬ 
ified. Namely, in the limit of k a <C 1—small solution 
salinities or large surface curvature 1/a- for the critical 
adsorption condition for PEs under confinement splits 
for the three fundamental geometries. We illustrate this 
behavior in Fig. [6] that is the main result of this study. 
Our results revealed a remarkable quantitative agreement 
with the recent theoretical predictions for the same sys¬ 
tems [37!]. The simulation approach enabled us to vary 
the polymer length and PE persistence, which are often 
quite problematic to be properly implemented from the 
first theoretical principles [18|. Also, we showed that for 
the critical adsorption onto concave surfaces more per¬ 
sistent chains require smaller surface charge densities to 
get adsorbed. For critical adsorption of PEs of varying 
polymerization degree N inside spherical cavities of ra¬ 
dius a, we found that the simulation results collapse onto 
a universal curve if the rescaled surface charge density is 
considered, namely a c -A a c aN. 

Finally, only the static properties of PEs under con¬ 
finement were considered in the current paper. It would 
be instructive as a next step to study the dynamics of 
charged polyions inside oppositely charged domains and 
cavities. In particular, the implications of polymer charge 
and adsorption to the spherical cavity interior can enrich 
the trends observed for looping kinetics of spherically- 
confined flexible and semi-flexible chains (87]. The video 
files of the Supplementary Material demonstrate, for in¬ 
stance, that the dynamics of PE chains in the adsorbed 
state is slowed down dramatically, as compared to de¬ 
sorbed configurations. This surface-mediated polymer 
confinement is consistent with the ultraslow relaxation 
of confined DNA molecules detected in single-molecule 
experiments during viral packaging (99j. Future devel¬ 
opments of the model will include the stud y of PE ad¬ 
sorption onto pH-responsive functionalized [1001 ] curved 
surfaces, the implications of a nonlinear ES potential dis¬ 
tribution on the position of the adsorption-desorption 
boundary a c (n,N) (93], and the adsorption of polymer 
chains with heterogeneous charge distribution. The lat¬ 
ter can be applied, for instance, to the surface-mediated 
adsorption of poly-peptide chains of partially folded pro¬ 
teins. 

Let us discuss some possible applications of our find¬ 
ings. Polymer encapsidation inside oppositely charged 


cavities m, m, m, m nnn is the fundamental mech¬ 
anism of assembly of cylindrical and spherical single- 
stranded RNA viruses [791, 0]. This process employs 
a delicately tuned adsorption of negatively charged RNA 
chains onto the positively charged interior of viral pro¬ 
tein shells. Direct applications of our observations to 
the properties of real RNA viruses might however re¬ 
quire the secondary RNA looped structure to be taken 
into account. The latter often pays an important role_in 
viral assembly and the RNA packaging process 0 . 1lO2| . 
Branching and self-association in the structure of com¬ 
pacted RNA yield, for instance, a weak overcharging of 
the entire virion: on average the negative charge of the 
enclosed nucleic acid chain is about 1.6 times larger than 
the positive charge of the enveloping protein shell [ 27 J . 

Another domain of possible applications includes the 
behavior of long DNA molec ules in micro-fluidic devices 
involving nano-channels |l03l . Il04| with attractive walls. 
Having in mind some appli catio ns to cylindrical channels 
of non-trivial cross-section (105| , one can consider in the 
future the PE adsorption on the interior of tubes with 
more complicated geometries, e.g. rectangular or trian¬ 
gular rather than circular cylindrical channels. Some ap¬ 
plications of our findings to the description of charge ef¬ 
fects of PE and DNA translocation though natural and 
synthetic nanopores Il06j | are also possible. One more 
immediate application of our results includes the prob¬ 
lems of protein adsorption—both in their native form 
and in the denatured state—in various porous media. 
For instance, polymer dynamics and adsorption in s ticky 
nano-channels of porous silicon studied in Ref. |l07l ] 
can pave the way for the selective separation of pro¬ 
teins from unknown mixtures, based on their surface 
charge and surface-adsorption properties. Moreover, de¬ 
fining the critical adsorption conditions is of vital im¬ 
portance for the fabrication of responsive and permeable 
multilayer capsids. They are being formed via the altern¬ 
ating adsorption of oppositely charged P Es [ 7 ] and used 
for diagnostic and therapeutic purposes [108| . Finally, 
nano-structured polymer-functionalized porous materials 
are used in electro-chemical super-capacitors em m3i 
and our results on polymers under extreme confinement 
might find some future applications in this area as well. 
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Appendix A: Supplementary Material 

In the Supplementary Material we include the video 
files illustrating the change of polymer conformations 
for the two basic geometries, as investigated in the 
main text. For each geometry, we fix the value of the 
surface charge density a and the confinement dimen¬ 
sions a and vary the solution salinity. In doing so, 
at small na values the polyelectrolyte chain is rather 
in the adsorbed state, while for larger na the electro¬ 
static polymer-surface screening gets stronger, the sys¬ 
tem crosses the adsorption-desorption boundary, and the 
chain desorbs from the interface (see also Fig. 6 of 
the main text). Parameters: the spherical cavity with 
the surface charge density of a = — 0.1/(47t) C/m 2 , the 
sphere radius is a = 50A, and the chain polymerization 
degree is N = 100, simulated at na = 1 (video 1) and 10 
(video 2). The cylindrical pore for the same values of a, 
N, and a simulated for na = 1 (video 3) and 10 (video 
4). Every video contains about 10 6 elementary simula¬ 
tion steps. Note that the length of the trace shown is 
about 0.1 of the entire trajectory length used in Fig. 6 
to determine the adsorption-desorption boundary. 
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